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and v[n] is the ECG signal. 

S9 Variation with Tpeak as mean evaluated in a symmetric interval of 10% of 
the Tstart-Tend-interval surrounding Tpeak, calculated by the formula: 

f Tpeak +0.Q5-(Tend -Tstart) ^2 

S9= £ (n - Tpeak) 2 • w[n] , 

^n=Tpeak-0.05(Tend -Tstart) J 

where w[n] = v[n]/ m 0 , 
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and v[n] is the ECG signal. 



S10 Variation with Tpeak as mean evaluated in a symmetric interval of 20% of 
the Tstart-Tend-interval surrounding Tpeak, calculated by the formula: 

f Tpeak+Q.\(Tend-Tstart) \ 2 

£10= Ya( n - T P eak f ' W W > 

\^n=Tpeak-0. \<(Tend -Tstart ) j 

15 where = v[/?]/ ra 0 , 

Tend 
n=Tstart 

and v[n] is the ECG signal. 
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S 1 1 The Hill parameter, K m , evaluated by least square fitting of the repolarisation 
integral, RI(t), from the Jpoint to the following Ponset as described by Kanters et al., 
"T wave morphology analysis distinguishes between KvLQTl and HERG mutations in 
long QT syndrome", Heart Rhythm (2004) 3, 285-292: 

RI(t) = K 



max 



S12 The Hill parameter, K m , evaluated by least square fitting of the repolarisation 
integral, RI(t), from Tstart to Tend analogous to the method described by Kanters et 
al., "T wave morphology analysis distinguishes between KvLQTl and HERG muta- 
tions in long QT syndrome", Heart Rhythm (2004) 3, 285-292: 



RI{t) = V t 
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Fl Flatness evaluated from Tstart to Tend, calculated by the formula: 

I 

f Tend \ 4 

F\= ^(n-m,) 4 -w[n] \ , 

\n=Tstart J 
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Tend - Tpeak 
where v[n] is the ECG signal. 



F12 Flatness parameter, Fl 1, normalized by the size of the R wave, calculated by 
the formula: 

F12 = 



v[Rpeak] - v[jpo int] | ' 



where v[n] is the ECG signal. 



F 1 3 Ratio of the total area under the T-wave from Tstart to Tend and the corre- 
sponding time interval, calculated by the formula: 

Tend 

ZvM 

JP^^ n=Tstart 

Tend — Tstart 
15 where v[n] is the ECG signal. 



F14 Flatness parameter, F13, normalized by the size of the R wave, calculated by 
the formula: 

FU 



20 F14 



v[/?/?ea£] — v[jpoint]| ' 



where v[n] is the ECG signal. 



25 F15 Ratio of the T wave height and the T wave width, calculated by the formula: 

v\Tpeak] 



F15 



Tend — Tstart 



where v[n] is the ECG signal. 



F16 The T wave height, calculated by the formula: 

F16 = v[Tpeak], 

where v[n] is the ECG signal. 



F17 Average slope from Tstart to Tpeak, calculated by the formula: 

F17 = v ^ T P eak ^~ v ^- Tstart ^ 
Tpeak — Tstart 
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A selection of these parameters is possible so that special genetic combinations are 
known with reference to the different stored parameters. The system can be updated by 
new data selected from different sources. 

5 The system and/or method can analyse the QT curvature of the ECG for indicating 
Long QT syndrome. This way, the Long QT syndrome can be indicated in an objective 
and effective manner which might occur in postsyncopal cardiac examination. 

The method can differentiate between different genotypes of the Long QT Syndrome, 
1 0 which is important for the treatment. It can, hereby, be achieved that the correct medi- 
cal treatments can be started. The system and the method can be used for test of drug 
influence on ECG curvature. 

The system can be trained, where the parameters' values are calculated for individual 
15 subjects, where an analysis of the parameters is performed such as a pattern classifica- 
tion method based on supervised learning, such as Discriminant Analysis, Nearest 
Neighbor Techniques, Multilayer Neural Networks, Decision Trees and Rule Based 
Methods or combinations of these. 

20 The final classification function is at least based on data from at least one LQT or drug 
influenced group and Normal subjects stored as a training set with the consequences 
that the classification method is improved by adding new subjects to the training set, 
which new subject can be tailored to demographic or gender differences. In addition it 
is achieved that reference values based on the training set can be selected from the 

25 most critical group of persons with reference to the parameters that are going to be 
tested. 

Once the parameters' values are calculated for individual subjects the mathematical 
analysis chooses the optimal (small) parameter set out of the complete set (large) from 
30 all categories, which values are stored as ref. values. It should be made clear that the 
final classification functions are based on data from at least one LQT or drug influ- 
enced group and Normal subjects (the training set) with the consequences that the dis- 
crimination method can be improved, in principle, by adding new subjects to the train- 
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repolarisation patterns. The most prevalent genes affected in LQTS patients are 
KvLQTl and HERG which account for more than 90% of LQTS genotype patients. 
The current study focuses on carriers of these two genes. Although some attempts 
have been made to develop quantitative measures that link different repolarisation 
5 abnormalities to specific LQTS related channel-opathities these methods have so far 
failed to provide a solid diagnostic yield. In current practice the duration of the QT 
interval is the only widely accepted quantifier of ventricular repolarisation. Yet, it has 
been recognized that the duration of the QT interval is only a gross estimate of repo- 
larisation since T-wave morphology is also important when characterizing the QT 

10 interval. This is evidenced by the fact that approximately 10% of all mutation carriers 
have a normal Bazett corrected QTc (<440ms) and 40% of KvLQTl and HERG carri- 
ers show QTc values between 410-470 ms that overlap with non-carriers. Conversely 
only 2% of all carriers present with a normal ST-T pattern and a normal QT interval. 
Morphological aberrations thus carry major implications for the identification of ab- 

1 5 normal repolarisation and have been included as diagnostic criteria equivalent to that 
of a positive family history for LQTS. 

Studies have shown that affected KvLQTl patients generally show broad based T- 
waves with a normal to relatively high amplitude and often without a distinct T - wave 
20 onset. For individuals with mutations involving the HERG gene the aforementioned 
studies have generally found low amplitude T-waves with bifid T-waves in 60% or 
more of the carriers. 

Cardiologists already include a qualitative assessment of T-wave morphology from the 
25 ECG in order to obtain information that augments the clinically established QT inter- 
val measurement and facilitates discrimination between LQTS genotypes. However 
qualitative description of repolarisation morphology may be biased due to intra- and 
interpersonal variability thus indicating the need for a standardized quantitative meas- 
ure of this parameter. 

30 

In the following is presented a novel multivariate categorization method that allows 
discrimination between KvLQTl, HERG and normal individuals based on Twave 
morphology recorded from 12-lead ECG's. Hallmark morphological features of T- 
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3. Results 

The discriminant functions were based on data from all KvLQTl, HERG and normal 
subjects. The 5 parameters included in both discriminant functions are listed in table 
2. 

5 The discriminative efficiency of both generated functions was statistically significant 
after inclusion of all 5 parameters (function 1: pO.OOOl, function 2: p<0.005). 

Variables Entered 



Step 


Entered 


1 


Fllstd 




QTcmeanVS 


3 


S5meanV5 


4 


D4std 


5 


S4meanV5 



10 Table 2. Variables used by the two discriminating functions. Stepwise introduction of 
more variables improved the ability of the functions to discriminate between KvLQTl, 
HERG and normal. 

A scatterplot was generated from the discrimination functions and groupings of indi- 
15 vidual genotypes can be seen in figure 4. The dotted lines were read from the SPSS 
generated territorial map and manually added. The lines reflect borderlines where the 
differences between each pair of discrimination functions are zero. All 16 processed 
ECG's were correctly classified and showed at least one discriminatory characteristic 
as defined by the 5 parameters included in the discrimination functions. Cross valida- 
20 tion of both discriminant functions was done with the leave-one-out method and all 16 
subjects were again correctly grouped. Reducing the number of variables resulted in 
misclassified cases due to lack of one or more discriminatory characteristics. In light 
of this finding we elected to perform further analysis of the selected parameters in 
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order to investigate the individual contributions of each variable to the separation of 
the three primary groups of subjects. Extreme values for all parameters were identified 
and the mean was computed. 

5 The result is plotted in figure 5. As expected the extent of interlead flatness variation 
observed in HERG and normal individuals was lower than that found in KvLQTl sub- 
jects. This is evidenced by the Fllstd parameter in figure 5a. When evaluating pa- 
rameter values S4meanV5 and S5meanV5 (figure 5b, d) the extent of asymmetry in 
KvLQTl and normal was generally less than that of HERG individuals. Both 

10 S4meanV5 and S5meanV5 are symmetry parameters and asymmetry in HERG indi- 
viduals was augmented in two ways: When bifid T-waves were present the interval 
from Tstart to Ttop was prolonged due to the definition of Ttop used in this study (the 
last highest point on the T-wave). Also, when the initial portion before Ttop was pro- 
longed in HERG individuals better discrimination was possible. Both phenomena 

15 were observed in HERG subjects. Generally the Bazett corrected QTc observed in 
HERG and KvLQTl was higher than that of normal individuals (figure 5e). However 
overlap existed between all three groups preventing separation of the groups by QTc. 
Since no single parameter included in the discrimination functions was able to sepa- 
rate KvLQTl, HERG and normal, we proceeded to investigate the classification effi- 

20 ciency provided by the three primary categories represented by the parameters in the 
functions. This was carried out by generating new discrimination functions using pa- 
rameters from one category only while excluding the other two. Then, from the new 
discrimination functions three additional functions were generated, this time allowing 
the inclusion of parameters from combinations of two categories. Scatterplots illustrat- 

25 ing the results of this analysis are shown in figures 6a-f. The first two functions (figure 
6a) included parameters that characterize the symmetrical properties of the Twave. 
83.1% of the 16 subjects were correctly classified. Arrows in figure 6a indicate the 3 
misclassified subjects. A second discriminant analysis was performed using flatness 
parameters. This resulted in 93.8% correctly classified subjects. Only one subject was 

30 not correctly classified as indicated by the arrow on figure 6b. The misclassified case 
was the same HERG subject incorrectly classified using symmetry parameters. The 
discriminatory efficiency of duration parameters was also evaluated. Discrimination 
analysis resulted in 93.8% correctly classified subjects. One HERG subject was mis- 
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Using only symmetry parameters, 3 subjects were misclassified. However no obvious 
visual characteristics on the three misclassified ECG's could be identified that ex- 
plained the incorrect classifications. The Bazett corrected QTc was 347ms for the 
normal subject, 425ms KvLQTl, 476ms HERG. Although an obviously prolonged 
5 QTc was present in the misclassified HERG subject it was not identified using sym- 
metry parameters alone. 

Discriminant analysis using parameters from the flatness category resulted in only 1 
misclassification. Again no visual characteristics were identified to account for the 
misclassification. Although it was anticipated that the 

10 

Figure 6. a) The result of discriminant analysis using symmetry parameters resulted in 
three misclassified cases (arrows). Visual inspection of the ECG's revealed no appar- 
ent abnormalities to indicate the reason for incorrect misclassification. b) The result of 
discriminant analysis using flatness parameters. One incorrectly classified HERG sub- 

15 ject was identified (arrow) even though no obvious visual abnormality indicated a dif- 
ferent genotype, c) Result of discriminant analysis using duration parameters. This 
result illustrates the failure of duration parameters to discriminate between KvLQTl, 
HERG and normal (arrow), d-e) Combinations of parameters from two categories 
illustrate the improvement in classification efficiency when compared to figures 6a-c 

20 evaluation of T-wave flatness would be able to discriminate HERG from KvLQTl 
subjects this was not accomplished by using flatness as a single descriptor of repolari- 
sation. Performing discriminant analysis based on the QTc parameter as the only vari- 
able resulted in 1 misclassification. This was not unexpected since it is well known 
that a substantial overlap in QTc values can exist between normal and affected indi- 

25 viduals. The lack of unambiguous discrimination between all groups by use of the 
QTc parameter alone emphasizes the hypothesis that additional parameters are needed 
to classify LQTS individuals. By combining parameters from two categories it was 
found that the discriminatory strength was increased. 

30 (figures 6d-f) This was evidenced by the fact that no subjects were misclassified using 
two categories. A particularly interesting finding, was the perfect separation of all sub- 
jects that was obtained using symmetry and flatness parameters with no duration pa- 
rameters included. This result implies the discriminatory strength inherent in parame- 



